Can the morphology of y-ray emission distinguish annihilating from decaying dark matter? 
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The recent results from the PAMELA, ATIC, FERMI and HESS experiments have focused attention on the 
possible existence of high energy cosmic ray e+e^ that may originate from dark matter (DM) annihilations 
or decays in the Milky Way. Here we examine the morphology of the y-ray emission after propagation of 
the electrons generated by both annihilating and decaying dark matter models. We focus on photon energies 
of 1 GeV, 10 GeV, 50 GeV (relevant for the FERMI satellite) and consider different propagation parameters. 
Our main conclusion is that distinguishing annihilating from decaying dark matter may only be possible if the 
propagation parameters correspond to the most optimistic diffusion models. In addition, we point to examples 
where morphology can lead to an erroneous interpretation of the source injection energy. 



I. INTRODUCTION 

Results from recent cosmic ray experiments (PAMELA fTl, 
ATIC|2|, FERMI L3J, HESSL4J) have raised the question 
of the origin of an "anomalous" population of high energy 
positrons in the Milky Way and motivated many studies. Cor- 
relation of the positron flux measured by PAMELA with the 
y-ray spectrum obtained by FERMI LAT is expected to give 
insight into the injection energy of the high energy electron 
and positron (e+,e^) population, and should also probe their 
spatial and energy distribution. In scenarios where high en- 
ergy e+,e^ are emitted by dark matter (DM), the spatial and 
energy distribution of this "additional" cosmic ray population 
is expected to follow the DM halo distribution at the injection 
energy E ~ Ei„j. This implies (assuming a spherical DM halo) 
that they should be spherically distributed with an energy den- 
sity that is maximal near the Galactic Centre. 
However, this picture could be modified if the high energy 
e+,e^ spatially propagate and lose energy in the galaxy ow- 
ing to inverse Compton and synchrotron losses. As a conse- 
quence of propagation, not only will the spatial and energy 
distributions of the high energy e+,e^ be modified but their 
final energy will be smaller than Einj. The y-ray spectrum 
obtained after propagation could therefore differ significantly 
from that obtained at injection. 

The issue of the y-ray flux associated with DM annihilations or 
decays into leptons has been addressed in several papers. For 
example, both the y-ray flux and y-ray spectrum in decaying 
and annihilating scenarios have been discussed in ref. ||5] |6l 
but propagation was actually neglected. More recently, the 
authors of ref.|7| have predicted the expected y-ray flux in a 
decaying DM model, taking into account e+,e^ propagation. 
Although computation of the flux is important, exploiting its 
value will be difficult owing to large uncertainties due to as- 
trophysical sources at these energies (Ref.|8|). Other papers 
have considered specific positions on the sky (e.g. interme- 



diate galactic latitudes, ||9l [TOl), or rely on very large-scale 
anisotropics (|7|). The work in ref. ifTTl raised the question 
of the morphology of the y-ray emission but mainly focused 
on the spectrum; however the propagation parameters adopted 
are not those favoured by MCMC studies (ifTHl. 
The question we raise in this Letter is whether the morphology 
of the y-ray emission alone (rather than the flux) can actually 
help to discriminate between the different DM scenarios. To 
address this issue, we compute y-ray maps originating from 
the interactions of e+e^ with the Interstellar Radiation Field 
(ISRF) spectra after propagation. We assume that the dark 
matter only annihilates or decays into e+e^ pairs and focus 
on y-ray energies that are accessible by the FERMI satellite, 
namely Ey= 1,10,50 GeV. We have neglected prompt y emis- 
sion which could arise from internal bremsstrahlung because 
its spectrum is model-dependent, however one should keep in 
mind that for some models, this emission could modify our 
conclusions. Given our assumptions, the injection energy of 
the e+ and corresponds to either the DM mass m^,,, or half 
the DM mass (£„,, — m^,„ or nidm/'^), depending on whether 
DM is annihilating or decaying respectively. We will consider 
three values of the injection energy: Ei„j = 100,500,1000 
GeV. At given Ei„j, the comparison between decaying and 
annihilating scenarios is immediate. In addition, since there 
are quite large uncertainties in the propagation parameters of 
cosmic rays, we will use three different sets of parameters re- 
ferred to as (MIN,MED,MAX) (cf Ref. 113 [H), which give 
a fair idea of the related uncertainty. We present difference 
maps of the y-ray contributions that highlight how morphol- 
ogy could help discriminate between the competing models. 



II. PRODUCING GAMMA-RAY MAPS 

To generate these maps, we apply the propagation scheme in- 
troduced by ifTTi , modified according to |fT5l[T6l for primary 
electrons. We compute the halo function / in terms of the 
electron energy: 
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TABLE I: Ellipticity for Ey = 10 GeV and Ei„j = 100 GeV for 0.1 ( 
subscript 10) and 0.5 ( subscript 2) of the intensity. 



where Qi„ = K /?„, is the Fourrier-Bessel transform of the 
source term (with K = ^ x (po/»idm)^ and = Ov the an- 
nihilation cross section if ^ = 2 and SI — F the decay rate if 
^ = 1; r| is the multiplicity, i.e. = 2 if Majorana particles, 1 
otherwise), L is diffusion slab half-thickness, Rgai is Galaxy 
radius and the propagation length expressed as: 

fEi„j 

Je„„„ b{E) 

where b{E) is the loss term and Ao,5 the diffusion parame- 
ters. Note that in Eq. ([T]i, r,z are cylindrical coordinates and 
a, are the roots of the Bessel function Jq. The gamma ray 
flux detected at the Earth is given by the integration along the 
line of sight of the convolution of the electron flux with the 
gamma emissivity of the electron interacting with the Inter- 
stellar Radiation Field (ISRF). The ISRF is mainly made of 
stellar light which is absorbed and re-emitted in the infrared 
by galactic dust. We use the model from 1 17J which can be 
fitted by a sum of blackbody (BB)-like spectra as in ITSl . Al- 
though this fit is valid only in the two kpc around the Sun, 
we expect that only the relative amplitudes associated with 
each BB vary from one position to another, but the tempera- 
tures should remain the same. The emissivity is computed in 
the same way as the losses in 1 18 1, making the approximation 
that the outgoing photon spectrum is a delta function ( lfT9l ) 
for each blackbody with which the electrons are interacting. 



m. RESULTS 

We now present the pixelized maps (with a pixel size of 1 
square degree) that we have obtained for the different sce- 
narios. Fig. [T] illustrates (for both annihilating and decaying 
dark matter models) the difference between the propagation 
patterns that arise by fixing Ei„j to 100 GeV and consider- 
ing three gamma ray energies Ey^l, 10,50 GeV. As one can 
see, in both cases, the e+,e^ which give rise to 1 GeV pho- 
tons have propagated further than those giving rise to 50 GeV 
photons. These features are common to all the maps includ- 
ing those obtained for heavier dark matter candidates. Note 
that the propagation parameters that we have considered to ob- 
tain this map correspond to the MED set (L = 4 kpc, 5 = 0.7, 
ifo = 0.0112 kpcVMyr). 



In Figj2j we have fixed Ey to 10 GeV and considered three 
values of (for both annihilating and decaying DM). In- 
terpreting the features for the particular case Ei„j = 500Gey 
and Ey = lOGeV is non-trivial. As can be seen from Table [H] 
because the ISRF is made of more than one BB, y-emission 
at 10 GeV can actually be due to more than one electron pop- 
ulation. Indeed bright emission at 10 GeV could be due ei- 
ther to electrons of ^20 GeV interacting with UV light or to 
^500 GeV electrons interacting with IR light. Hence as seen 
in Fig|2] the 10 GeV emission is nearly spherical, and could 
be interpreted either as an injection energy of ^20 GeV or of 
^500 GeV, leading to very different interpretations concern- 
ing the mass of the DM particle. However this degeneracy can 
be lifted by looking at higher energies, as electrons injected at 
20 GeV cannot produce gamma rays of 50 GeV. This thresh- 
old effect stresses how important it is to look at different y-ray 
energies and to compare the various morphologies in order to 
understand the properties of the DM. 

In the third column, we exhibit the difference between the two 
normalized maps (decays — annihilations) so as to exhibit the 
differences of morphology between these two emission mod- 
els. The negative values at the Galactic Centre confirm that 
the e+,e^ from annihilating DM are mainly produced locally 
and that propagation cannot completely smooth out the con- 
trast with respect to decaying DM electrons and positrons. 
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TABLE II: Electron energy E^ responsible for the emission of an 
electron of energy Ey through inverse Compton scattering on each 
blackbody component of the ISRF. 



In Fig|3] we show the effect of the propagation parameters for 
Ei„j ~ 1 TeV. As is expected, the e+e^ diffuse far more for the 
set of propagation parameters MAX (for which L = 15 kpc) 
than for MIN. Although it may be possible to constrain de- 
caying versus annihilating DM in the MAX and MED cases, 
it seems impossible to distinguish these two scenarios in the 
MIN case. To compare the propagation features between an- 
nihilating and decaying scenarios, it is useful to look at the 
ellipticities Eio and 82 of the y-emission. To define these quan- 
tities, we measure the size of the semi-major axis oio (or 02) 
and the semi-minor axis biQ (or b2) of the ellipse that has an 
intensity of one tenth (or one half) of the maximal intensity. 
Ellipticity is then defined as 1 — b/a. The results are summa- 
rized in the caption of Fig. l|2|3| and Table [ll] Only for larger 
masses, or in the optimistic case where the sensitivity allows 
us to measure ejo, is discrimination possible, especially for 
the MAX propagation model. 
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FIG. 1: Annihilating versus decaying DM for Ei„j =100 GeV and Ej= i, 10, 50 GeV. In tliese figures, fluxes are normalized to the central bin 
so as to make the comparison of propagation length obvious. EUipticities at 0.1 of the central bin intensity are equal to eio = 0.41,0.31,0.15 
and eio = 0.75,0.72,0.68 for annihilating versus decaying DM respectively. 




FIG. 2: Annihilating versus decaying DM for Ej = 10 GeV and Ei„j = 100,500, 1000 GeV. EUipticities at 0.1 of the central bin intensity are 
equal to Eiq — 0.31,0.06,0.35 and Eiq = 0.72,0.67,0.72 for annihilating versus decaying DM respectively. 



IV. CONCLUSION 

We have generated maps of y-ray emission associated with 
e+,e^ population originating from DM annihilations or de- 
cays. We show that propagation is important for both DM 
scenarios, but although the propagation features differ, they 
are difficult to distinguish if the propagation parameters corre- 
spond to MIN (and perhaps MED) rather than to MAX. This 



is, in fact, surprising, as one might have expected these two 
scenarios, which involve distinct powers of the dark matter 
density, to differ significantly. Actually, in the MIN case, de- 
tection would be extremely challenging since most of the sig- 
nal would be hidden by galactic sources. In some cases, the 
IRSF can make the Galactic Centre bright enough to be mis- 
interpreted as e+e^ with a lower injection energy. 
We have verified that changing the energy density of the ISRF 
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FIG. 3: Annihilating versus decaying DM for Ei^j =1 TeV and the min and MAX propagation parameters. Ellipticities at 0. 1 of the central bin 
intensity for MIN and MAX are equal to eio = 0.67, 0.35, 0. 15 and eio = 0.92, 0.72, 0.58 for annihilating versus decaying DM respectively. 



has little effect as the increase of the y-ray emissivity is par- 
tially compensated by the electron density decrease due to in- 
creased energy losses. Varying the intensity of the magnetic 
field within reasonable values has also little impact as syn- 
chrotron emission is not the main energy loss term in most 
cases. In both cases the impact is mainly on the intensity and 
not on the ellipticity. However a full spatial description of 



both the ISRF and the magnetic field could have effects that 
are beyond the scope of our analytical approach. 
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